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SUMMARY 


A simplified analysis is made of mass transfer cooling - that is, 
injection of a foreign gas - near the stagnation point for two-dimensional 
and axisymmetric todies. The reduction in heat transfer is given in 
terms of the properties of the coolant gaB and it is shown that the heat 
transfer may he reduced considerably by the introduction of a gas having 
appropriate thermal and diffusive properties. The mechanism by which 
heat transfer is reduced is discussed. 


INTRODUCTION 


The reduction of heat transfer near the stagnation point of a blunt 
body (two-dimensional or axisymmetric) is of primary importance when the 
body attains high velocity. One method of cooling iB to introduce gas 
of high specific heat Into the laminar boundary layer near the stagna- 
tion point and to allow this gas to flow over the nose of the body so 
that large amounts of heat are convected away from the nose. 

A number of exact solutions based on various assumptions of fluid 
properties are available for stagnation-point heat transfer with no 
injection. (A comprehensive list of these references may be found in 
refs. 1, 2, and 3*) It has also been shown that the variation of the 
product of density and coefficient of viscosity across the boundary 
layer has an important effect on the heat transfer (ref. 4 ). 

The boundary -layer equations for a binary mixture are well estab- 
lished (ref. 5 ) and- several methods are available for the evaluation of 
the thermal and diffusive properties of a binary mixture (ref. 6 ) . 

Exact solutions (refB. 1 and 3 ) have recently been obtained for 
air-to-air injection near the stagnation point; these solutions have 
required lengthy numerical integration techniques and are necessarily 
limited to discrete values of the parameters involved (for example, the 
rate of injection of mass and the viscosity-law assumption) . 
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The purpose of this paper is to present a simplified analysis "by 
which the effect on the heat transfer of injection of air, or a foreign 
gas, may be found and, also, to show more clearly the important thermal 
and diffusive properties of the foreign gas if it is to be effective as 
a coolant. Application of the present method to ablation cooling is 
made in reference 7 * 


SYMBOLS 


x 

y 

z 

Y 
Z 
u 

V 


V 

U 

C 

T 

T 

W 

W 

P 

d 

k 


coordinate along wall 
coordinate normal to wall 
transformed y-coordinate 

arbitrary value of y outside boundary layer 
arbitrary value of z outside boundary layer 
component of velocity in x-direction 
component of velocity in y-direction 

modified velocity component, — 

^w 

free -stream velocity in x-direction 
constant in velocity distribution 
temperature 

mean temperature (eq. (43)) 
concentration of foreign gas 
mean concentration (eq. (42)) 
density of mixture 
coefficient of viscosity 
thermal conductivity of mixture 
coefficient of binary diffusivity 
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m 


°P 

°P 


Hi 


Nu 


H- 


Pr 


N, 


Sc 




rate of mass injection per unit area of wall 
specific heat at constant pressure 

mean specific heat (eq. (k'j)) 

Nusselt number 

Reynolds number 
Prandtl number 

Schmidt number 

heat -transfer rate per unit area 
velocity-boundary-layer thickness 


8 t thermal -boundary-layer thickness 

by concentration-boundary-layer thickness 

K function of Prandtl number and Schmidt number (defined by 

eg.. (10)) 

Heff effective heat capacity of foreign gas 

Subscripts: 

e external flow, near stagnation point 

w wall 

c coolant, far inside wall 

0 no injection 

1 foreign gas 

2 air 
b 


body 
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DESCRIPTION OF THE BOUNDARY LAYER 


The flow considered is that shown in figure 1. In the steady state 
the laminar boundary layer near the stagnation point is a thin layer of 
fluid in which the velocity, temperature, and concentration of the 
foreign gas vary rapidly from the external stream values to the wall 
values. In the neighborhood of the stagnation point the three super- 
imposed boundary layers (that is, the velocity, temperature, and con- 
centration boundary layers) have constant, although different, 
thicknesses. 

It is assumed that the component of velocity parallel to the wall 
is linear in x and that the normal component of velocity, the tempera- 
ture, the concentration of the foreign gas, and the properties of the 
mixture are all functions only of y, the distance normal to the wall. 
(This assumption results directly from the "similarity" nature of the 
flow and the absence of thermal and concentration gradients along the 
wall.) It is also assumed, consistent with the foregoing variation of 
velocity components, that the coolant gas is injected normally at the 
wall with a velocity independent of x. 

When a given amount of coolant gas is injected, it diffuses through 
the boundary layer and is convected with the air, as a mixture, under 
the action of the pressure gradient imposed by the external flow and 
the shearing stress due to the presence of the wall. The concentration 
of foreign gas at the wall is uniquely determined by the rate of injec- 
tion. The extent of the boundary-layer shielding depends upon the spe- 
cific heat of the foreign gas and upon the wall temperature as well as 
the coolant and stagnation temperatures. The shielding also depends 
upon the manner In which the coolant gas diffuses through the boundary 
layer. The analysis shows the relative importance of these quantities. 


AMLYSIS 


The method to be used is as follows: the wall conditions for 

coolant injection are formulated In a simple manner (and the results 
are justified by a more detailed consideration of the flow within the 
porous wall presented in the appendix) ; use is then made of previous 
"exact" results for stagnation-point heat transfer with no injection in 
order to determine the gross characteristics of the thermal and viscous 
boundary layers near the stagnation point. These results are then used 
In a simple approximate integral method to determine the effect of injec- 
tion on the heat transfer to the wall. 



MCA TN 4391 


5 


The Coolant Flov 

The "boundary conditions at a porous wall through which a coolant 
is injected may he obtained very simply by ignoring the presence of the 
solid part of the wall; the justification for this is given in the 
appendix where more detailed considerations are made. 

When the steady flow of coolant toward the surface of the wall is 
considered, the volume taken up by the solid wall being neglected, it 
is seen that there is a balance of diffusion and convection which 
governs the flow of mass and heat within the wall. 

The transfer of mass is given simply by 

prv = p^Vy = m (l) 


where m is constant. 

The diffusion of air inward from the surface is balanced by the 
convection toward the surface: 

-P D 12 || = d - W)m (2) 

Diffusion of air Convection of air 

from surface toward surface 

since W is the concentration of the foreign gas and (l - W) is that 
of air. 

Similarly, the transfer of heat is given by 

k ^ = c p,ld " T c)“ (3) 

Diffusion of Convection of heat 

heat from toward surface 

surface 

It is important to note here that, even though the specific heat of the 
mixture Cp is given by 

c p = c P ,i w + c P , 2 d - w ) . d) 

(with Cp 2 . and- Cp 2 constant), the value used correctly in equation ( 3 ) 
is c p,lJ unit time the heat transferred at any point y in a 
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direction away from the surface is sufficient to raise an amount of 
foreign gas of specific heat c p,i through .the temperature range T - T c . 

Evaluating equations (2) and (3) at the wall surface gives the fol- 
lowing boundary conditions which are required for the solution of the 
boundary-layer equations: 



= (1 - W w )m 

(5) 

w 


p,l( T W ** T c)® 

(6) 


Detailed study of the flow within the porous wall yields the same bound- 
ary conditions. (See appendix.) 


The Boundary -Layer Integral Equations 


The integral equations which describe the transfer of mass and heat 
in the boundary layer are derived in a simple way without reference to 
the general differential equations. 

Consider a small rectangle of height Y (Y greater than any 
boundary -layer thickness) and length Ax, in which p, W, and T are 
independent of x, near the stagnation point of a two-dimensional body. 
(See fig. 2.) The continuity of mass may be expressed aB follows; 


-p(Y)v(Y) Ax + m Ax 



pu dy 


Jx+Ax 


Flow in. from Injection 

external stream 


Flow out at 
x + Ax 


This equation, in the limit as Ax -» 0, becomes 


r Y 

-p(Y) v(Y) + m = A. / pu dy 

tix J 0 



pu dy 


Jx 


Flow in 
at x 


(7) 


When the form TJ = Cx of the external velocity is used and it is noted 
that 
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and 


u = C - x 
U 


where ^ is independent of x, then equation (j) reduces to 


-p(Y)v(Y) + m 




ay 


( 8 ) 


Since the density p is a function of concentration and temperature, 
it is convenient to introduce the quantities 


pv 

V = ^ 


dz = — dy 


( 9 ) 


so that equation (8) becomes finally 


r z 

-P^VCZ) + m = Cp w J g dz 


( 10 ) 


Similarly, the equation for the transfer of foreign gas becomes 


m 


0p "/o 


dz 


( 11 ) 


Gas injected 


Gas convected 


The equation of heat transfer is written 
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-c p>2 (T e - T W )P W V(Z) 


Convection of heat 
from external stream 



Transfer of 
heat to wall 



( 12 ) 

Convection of heat 
in boundary layer 


The quantity V(Z) may he eliminated from equation (12) by use of 
equations (4), (10), and (ll) j thus, 

CP W / Q Z Pp(Te - T) 2 a. . A + c pA (Te - T„)m (15) 

Alternatively, by use of the boundary condition given by equation (6), 
equation ( 13 ) has the form 



T) g dz = c p ^ L (T e - T c )m 


(14) 


The equations for transfer of heat and mass near the axisymmetric 
stagnation point can be derived in a similar manner and are identical 
with equations (10) to (l4) where x is again the distance from the 


stagnation point along the surface but 


z 



dy. 


Heat-Transfer Relations for No Injection 

The rate of heat transfer is first given in terms of Nusselt num- 
ber Reynolds number and Prandtl number Np r ^ v for which 

the usual definitions near the stagnation point are used: 


%u,w 



R w = 


tyjx 

^w 


X 2 
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and 


l^f 

N Pr,v = c p,2 


The rate of heat transfer to the vail with no injection is now 




(* fl ’ ^ - M^) 172 


(15) 


r Ni 


Exact evaluation of the dimensionless parameter 


Nu^v' 

1/2 

) 0 


has been 


made hy many investigators (for example , refs. 1, 2, and 4). The 
following formulas were found to he in good agreement with previous 
exact results (ref. l) under the assumption of constant pp: 

For the axisymmetric body. 


N- 


Nu,w \ 

A 1/2 jo = 


0.765 - 0.065 [l - 


N- 


Pr,v 


0.4 


(16) 


and for the two-dimensional body, 


3 


Nu,w \ 

1/2 

0 *- 


0.570 - 0 . 065(1 - ^ 


N- 


Pr,w 


0.4 


(17) 


Alternative expressions, based on the solutions presented in refer- 
ence 4 in which pp is not constant, are as follows: 

For the axisymmetric body, 


= 0 . 765 ^ 


>0.4 
I N- 


0.4 




Pr,w 


( 18 ) 
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and for the two-dimensional body. 


^ 1/2 /o 


=0.570/^ 


tPr,w 


Boundary -Layer Thickness Without Injection 

Use is now made of these results (eqs. (16) to (19)) to determine 
the thickness 8 U q of the velocity boundary layer without injection. 

Equation (13) with m = 0 is written in the form 

1 r ZT e- T »,. _ 1 AM 1 / 2 1 /V w\ 


1 r Z T e ~ T u dz _ 1 /IV\ X / 2 1 |^Nu,w\ 

u,0 '-I 0 T e “ T w u s u, 0 \^Cy %r,wL, l/2 1 


and linear profiles are assumed for velocity and temperature in order to 
evaluate the integral, in equation (20); that is. 


U 5,. 


T e - T 
T e ~ T v 


= 1 - 


®T,0 


where 8^, 0 is ' thiGlm ess of the thermal boundary layer. The left- 

’ , , ifac) 2 

hand side of equation (20) is then approximately — L - - and equa- 

6\pu,oy 

tion (20) reduces to 


l/ S T,o\ 2 = 1 f^) 1 ^ 1 

6 \&u,C / s u,0\ p w C / %r,wL 1 / 2 
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If the value 
script (l). 


N- 


of 


Nu,w \ 

> l/2 jo 


when = 1 is signified by the super- 


then equations (l6) to ( 19 ) show that 




(1) 

0 


%r,w 


OA 


(23) 


This equation shows the effect 
coefficient. When' Npj. w = 1, 


are similar so that 


^,0 

8 u,0 


= 1; 


of the Prandtl number oh the heat-transfer 
the thermal and viscous diffusive effects 

thus, from equation (22) 



(2lt-) 


The nondimens ional thickness 


^,0 



is now assumed to be inde- 


pendent of the Prandtl number (since the momentum equation is only 
weakly dependent on the Prandtl number); thus, when equations (23) and 
(24) are inserted into equation (22), the following simple result is 
obtained: 


S T,0 

8 u,0 


= N- 


Pr,w 


-0.3 


This result is analogous to that for a flat plate (ref. 8): 


( 25 ) 


s T 

— = N- 
8 u 


'Pr,w 


-1/5 
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Injection of Air 

Consider now the effect of injection of air into the "boundary layer. 
The "boundary layer "becomes thicker because of the increase in mass flow 
which results in reduced gradients in the boundary layer. In particular, 
the heat transfer to the wall is reduced since the heat convected par- 
allel to the wall in the boundary layer is increased. 

The effect of air injection on the boundary- layer thickness is 
given by considering the mass flow in unit time per unit area, an 
follows : 




+ 


m 


Mass flow with 
injection 


MaBS flow with Mass 

no injection injected 


which gives 




m 



( 26 ) 


when linear profiles are assumed for ^ ^that is, ^ 
tion) and — = -5- (injection)) . 

u &u 7 


— — (no injec- 
5 u ,0 


It is assumed that the relation given by equation (25) is still true 
since the ratio of thermal-boundary-layer "thickness to velocity-boundary- 
layer thickness depends on diffusive processes, that is, on the Frandtl 
number . 


With Cp ; i = °p,2 equation ( 13 ) may be written by use of 
equation ( 15 ) In the form 


f Zl e - T u d2 . if.) 

1/2 

1 %u,w m 

... - -J 

J 0 T e - T w U 5 u\Pw C y 
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Linear profiles for temperature and velocity (eqs . (21)) and the expres- 
sions for 6 U (eq. ( 26 )) and 5 Uj0 ( ec l‘ ( 2l{ -)) are use<i to simplify 

this equation, with the following result: 


1 N Nu,w _ 1 /^Nu,w\ 

^ c 75 = ^wi 0 


1 - r N- 


- « -°A 

3 P^w j 


m 


(fW 3 ) 


1/2 


(27) 


In terms of the rate of heat transfer q, equation ( 27 ) has the 

form 


<1 


*0 


Heat transfer 
for no 
injection 



°S,2&e 



Boundary-layer shielding 
by convection 


The quantity ^1 - i Np r ^ w -0 ‘^j (T e - T w ) may be interpreted as the 

average temperature rise in the boundary layer of the mass introduced 
at the wall surface. An alternative expression of this heat balance, 
obtained by using the boundary condition given by equation (6), is 


% 


= c 


P,2 


(T w - T c ) + 



-0.6) 


(Te - Tv) 


m 


(29) 


When q^ is expressed in terms of dimensionless quantities by use 
of equations ( 15 ) and ( 23 ), equation (29) becomes 


T 

T 


e 


e 



N- 


Pr,w 


0.6 


m 


(PwPw C ) 


1/2 


N: 


Nu,w\ 

m 


( 1 ) 


m 




5 (p»t%c) 1/2 


( 30 ) 
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Equation ( 30 ) may be used to determine the rate of injection 
required to maintain the wall temperature T w at a given desired value 

when T e and T c are specified. 


Before proceeding to considerations of injection of foreign gas, 
the results given by equation ( 27 ) are compared with available exact 
solutions as a check on the validity of the several simplifying assump- 
tions that have been made. 


Firstly, under the assumption of constant pp, equations (l6) 
(eq. ( 17 ) for the two-dimensional flow) and ( 27 ) are used to give 

Njju . 

— as a function of the dimensionless rate of mass injection 

«w 


m 


(Pw^C) 


1/2 ‘ 


A comparison of these results with the results of refer- 


ences 1 and 3 shows very good agreement except for the extreme rates 
of injection with Prandtl number equal to unity. (See figs. 5> 4-, 
and 5 .) 


Secondly, for variable pp, equations (18) (or eq. (19)) and ( 27 ) 
are used and the results are compared with those of references 3 and 4. 
For a Prandtl number of 0.74 and no injection reference 4 gives, for 
the axisymmetric case. 



0.67, 




e^e' 




0.4 


and the present report 


gives (from eq. (l8)) 


%hi,w 



O .667 



A comparison of the results obtained by use of equation (18) with the 
results of reference 3 in which the Sutherland viscosity law was used 
is presented in figure 6. 

The effects of injection are also compared when pp is variable; 
again good agreement is found except for the high injection rates . 

(See fig. 5 .) 
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With this justification of the assumptions that have been made, the 
method is now extended to take account of the additional effects which 
result from the injection of a foreign gas (that is, gas having prop- 
erties different from those of air) . 


Injection of Foreign GaB 


The mechanism by which air diffuses through the boundary layer in 
a binary mixture due to the concentration gradient is exactly analogous 
to that by which heat diffuses In the absence of viscous dissipation 
for air-to-air injection. In this analogy the coefficient of binary 

diffusivity Dip corresponds to the thermal dlffusivity and 

PC P ,2 


the Schmidt number 



PC P ,2 

k 


Ng c = - y — 
3 PDi 2 


corresponds to the Prandtl number 


The analogy is easily seen when the gas-transfer equation (ll) is 
rewritten as 


cpu T Z JL H dz = A 
^ Jo W w u W w 


(31) 


and compared with the heat-transfer equation for air-to-air injection. 


1 Z r p _ m _ 

, 1 e x u , m 

Cp„ / — dz = 


vJ 0 


T e - T w U T e - ^ 

T e “ T c 


( 32 ) 


The boundary condition at the wall Is 


- PD12 = U “ W w )m 


( 33 ) 


which is compared with 


k dT 


pc p ,2 dy 


= (T v - T c )m 


( 34 ) 
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The analogous quantities are presented in the following table: 


Binary 

diffusion 

Thermal 

diffusion 

W 

EH 

1 

0) 

EH 

T e ~ T c 


T - T_ 

1 - w 

c 

T e " T c 

Ww 

T e - Tv 

T e - T c 

D 1 2 

k- 

pC P,2 

N Sc 

N Pr 


It may he verified that all the boundary conditions at the wall, in the 
external flow, and at large distances into the interior of the wall 
agree with the foregoing analogy. 


Because of the foregoing comparison the following results may be 
written immediately: from equation ( 25 ), 


= N -0.3 

5 U Sc,w 


(35) 


where is the thickness of the concentration boundary layer, and 

from equation ( 30 ) 


N, 


Sc,w 


0.6 m 


W w “ 


(p^C ) 1 / 2 


'W\ (l) 1 


m 




1/2 


(36) 


0 
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The value of 


/ %u,y\ ^ 

U v 7o 


to be used in equation ( 36 ) is necessarily 


that given by equation (18) or equation ( 19 ) since p|j. is now a vari- 
able quantity and depends on the concentration of foreign gas in the 
boundary layer. 


When the gas injected is different from air, the amount of shielding 
by convection in the boundary layer depends critically upon the specific 
heat of the gas c . 

The heat-transfer equation ( 13 ) is rewritten, by use of equations (4) 
and (6), as 

CfV / 0 T° P ’ 2 + fa’ 1 ‘ °P,2) W ] (Te ' dZ “ c PA 4(Te - T = ) (37) 

and the linear profile 


W 


= 1 



(38) 


is assumed. 


Substitution of equations (21), equation ( 38 ), and the definition 
of % u into equation ( 37 ) gives results similar in form to those for 

air-to-air injection (eq. ( 27 )): 


N n 


N. 


Nu,w _ 1 / Nu,w 

Pr,w b // 2 N Pr,w\R v 1 / 2 / 



(PwPw^) 


C PA _ 

°P,2 


11 - KN, 


Sc,w 


0 . 6 \ 


m 


(Pw^ c ) 


172 


(3 9 ) 
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where 


K = 6 — f Z JL 
s u J 0 W w 


T e - T u 
— £ dz 

T e " U 


(40) 


and is a function of Np r>w and %c ; w' 

Equation ( 39 ) may he expressed in the alternative form 


* - «o - 




KNsc.w 


0.6 


1 - i %r , W ‘°- 6 


°P,2 


KNsc, 


0.6 


1 Nt.- 0 - 6 
5 H Pr__. 


1 “ 5 ^Pr ,w 


- 0.6 


(* “ 5 N Pr,v'°- 6 )( T e - T V )A (41) 


The second term on the right-hand side of equation (4l) has the form 

[°p,lM + «p, 2 U- - «>] - Tvlm 

where W, the mean concentration of the foreign gas in the boundary 
layer, is 


if . 1 • g sc,v 
1 - 5 W 


0.6 

■ 0.6 


(42) 


and where T, the mean temperature in the boundary layer, is 


T - T w + (l - \ N Pr,w"°' 6 )( T e - T w ) (43) 

In terms of these mean values equation (4l) becomes 

q = qo [=p,lW + c p ^ 2 (l - W)] (T - T w )m 

( 44 ) 

Heat Heat transfer Boundary-layer shielding 

transfer for no 

injection 
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Thus, in order to achieve the maximum amount of shielding for given 
temperature conditions and mass injection, c p^l should he made as 

large as possible and, in addition, T should he made as large as possi- 
ble; if c_. 1 > c rs then W should he large, and if c_ < c 0 


then W 

£ - ?W 

T e - 

that W 


should he small. The variations with N pr ^ w and N gc ^ w of 
and W are shown in figures 7 and 8, respectively. It is seen 

T - T„ 

increases as „ decreases and that — increases as 

bc,w m _ m 


Npr v increases . 


Maximum shielding in the boundary layer is therefore achieved when, 
for c p -l > c p 2 , N Pr ,w is large and N Sc, w small and, for 

c p,l < c p,2 J N Pr,w ls ^se and Ng^ is large. 


These diffusive effects on the shielding are explained qualitatively 
in the following way (fig. 9 ): the convective shielding is most effi- 

cient when the gas of higher specific heat is transported in the regions 
of highest temperature and velocity, that is, in the part of the boundary 
layer farthest from the wall. 

Equation ( 25 ) shows that when is large the velocity boundary 

layer is thicker than the thermal boundary layer and thus more of the 
hot gas mixture is convected. When > c -p,2 -^Sc,w is sma ~ 1 ~ l > 

equation ( 35 ) shows that the concentration-boundary-layer thickness is 
greater than the velocity-boundary-layer thickness - that 1 b, the foreign 
gas of higher specific heat diffuses quickly through the boundary layer 
before being convected. On the other hand, when Cp,i < Cp^ g an< ^ ^Sc,w 
is large, the foreign gas remains near the wall and displaces the air of 
higher specific heat into the hot fast region where it is convected. 

An important consideration in any cooling system is the weight 

penalty which is incurred in order to maintain the wall at a desirable 

temperature. The relation between the rate of mass flow m and the 

temperatures T e , T w , and T c is obtained by substitution of the 

expressions 


0 . = Cp,i(T w - T c )i 
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and 


90 - =p,2 (I « - Tw)t(W:) l/2 

(from eqs. (6) and ( 13) ^ respectively) into equation (4l) . This suhsti- 
tution results in the following expressions f^vhen either equation (18) 


or equation (19) is used for f — ^ u . ; v \ 


*v l/2 J 0) 


For axisymmetric flow. 


°* 765 (p^) N Pr,w °‘WwC) ^ c p ^ 2 (T e - T w ) 
m = — L „ — 


®p,l(^-v - T c ) + Cp(T “ T w ) 


(45) 


and for two-dimensional flow, 

, 0. 14- 


0.570 


'P e ^e' 

iPw+Hr; 


%r,w (PwPw^) ^ ^,2^6 


- T w) 


m = 


c p,l(^w " T 0 ) + 5p(T - T w ) 


m 


where 


5 P - °p,l w + c p,2( 1 “ W) (47) 

Prohahly the most important parameter is the total effective heat 
capacity of the coolant; this parameter is a measure of the total 

heat absorbed when unit mass of the coolant is used and is defined as 
H eff “ y “ c p,l( T w - T c ) + Cp(T - T v ) (48) 
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METHOD OF CALCULATION OF RATE OF MASS INJECTION 


The problem of most interest from practical, considerations is that 
of determining the rate of mass injection m necessary to achieve a 
desired wall temperature T v when the stagnation temperature T p and 

the coolant temperature T c are given. This calculation is made com- 
plicated "by the dependence of the mean specific heat c p (required for 

eq. (45) or eq. (46)) on the rate of mass injection; both these quantities 
depend on the concentration W v . (See eqs. (4-7) , (42), and (35)-) 

The following method, however, gives results with relatively little 
computation: 

(1) The parameters Np^, w and N gc v are found in terms of W w 

for the particular binary mixture under consideration. (Simple empirical- 
methods are given in ref. 6.) 

( 2 ) Equation (43) 


T ^ = x 

T e - T w 


1 

3 



6 


Is used to give T (fig. 7) • 

( 3 ) Equation (42) 


W = 


1 - KNsc,w 


0.6 


1 - I %r,v' 


- 0.6 


is used to give W (fig. 8). 

(4) Equation (47) 


% = °p,l W + c p,2^ ' 


then gives c p (fig. 10) . 
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( 5 ) A plot of H rVM against W w can now be made 

(p^C^Wej 

from the results of the substitution of equations (43) , (42), and (47) 
into either equation ( 45 ) or equation (46) . 


( 6 ) A second plot of 2 — — 7 -I 1 against W w can be 

(p w ^i w C) 1 / 2 \ p e^e/ 

obtained as a result of converting equation ( 36 ) into the following form 
for the axisyrametric case: 


m 


(fV^ c ) 


/PyPy^ 

1 / 2 ve^e , 


0.4 


0.765tfw 


N, 


Sc,w 


0.6 1 y 

“ J "W 


(49) 


or into the following form for the two-dimensional case; 


m 


/Pw Py 


0.4 


O. 570W w 


(<W C ) 


l/2 v e P e> 


N 0.6 1 y 

N Sc,v ” 5 w w 


(50) 


The intersection of the curves found in steps (5) and ( 6 ) gives W w and 
0.4 


m 


(PyPyO) 


fPyPy ' 


•^/2 \PePe ) 


( 7 ) Now, Py-tV can be found from W w and T w , and p e |i e from 
external conditionsj thus, m can be calculated. 

Step (7) may be omitted if the approximation 


' p yPy \ 


0.4 


' p yPy 

\ p e p e/ 


1/2 


is made, in which case 


m 


( PePe^) 


1/2 


is given directly by steps (5) 


and (6) . 
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In order to see the main effect of 


■» 

”Pj2 


on the mass -flow require- 


T - T 

ments for given temperature parameter — values of the Prandtl 

T e - T v 

number and Schmidt number equal to unity are taken. 

Thus, from equation (43) 

* - 2 
T e - T v " 3 

from equation (40) K = i, and from equation (42) W = — . Using this 

2 4 

value for W in equation (47) gives 

Cn = c_ , + i c „ 

P 4 P ,1 4 P,2 

Then, equation (45) for axisymmetric flow becomes 


m 


m 


(p e P e c ) 


1/2 


( p w tJ v ! - : ) 


X / 2 \pePe / 


0.4 


0.765 


T C + l\ + 1 


®P,2\T e - T w 2/ 6 

and equation (46) for two-dimensional flow becomes 


m 


0.570 


( PePe 1 -") 


1/2 <=p,nA, - T c x 


+ - 1 + i 


^,2 N^e - 2 J 6 


The results are shown in figure 11. It is seen that the required 

c 

mass flow is reduced by a large factor for the higher values of fij x . 

°P>2 

This is especially true for the higher heating rates fas — — -» 0 ( 

V T e - ^w ) 

inasmuch as the shielding by heat convection in the boundary layer is 
large . 
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CONCLUDING REMARKS 


An approximate analysis has been presented whereby the reduction 
in heat transfer near the stagnation point may be calculated with little 
difficulty when the coolant properties are known. The agreement with 
available exact solutions for air-to-air injection is extremely good 
in view of the simple approach employed, and the qualitative trendB of 
the results for foreign-gas injection in explaining the shielding mecha- 
nism suggest that the approximate analysis will generally give reliable 
results. It is expected that the simplified calculation which shows 
the dependence of the coolant mass requirement on wall temperature and 
coolant specific heat will provide a good estimate for engineering 
purposes . 

The conclusions of the analysis are summarized briefly as follows; 

Maximum boundary-layer shielding is achieved when the gas of higher 
specific heat is convected in the hot fast-moving part of the boundary 
layer farthest from the wall; this requires that the coolant gas intro- 
duced have 

(a) High specific heat compared with that of air 

(b) High binary diffusivity (that is, small Schmidt number) 

(c) Large Prandtl number . 

Of these three properties the first is the most important, whereas 
the diffusive effects are of secondary importance. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 3, 1958. 
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APPENDIX 

FLOW WITHIN TEE POROUS WALL 


For the purpose of the present analysis it is assumed that the wall 
is uniformly porous and presents , to the stream, an area A of solid 
material and an area B across which coolant gas is transported. 

The manne r in which gas is transported through the porous material 
is given by 


-p Dl2 g B 

Diffusion 
toward surface 


-t- pvWB 

Convection 

toward 

surface 


m(A + B) 

Total gas 
transported 


(Al) 


since m is the rate of mass introduced per unit area of body surface. 

The air contained in the mixture is transported according to the 
following equation: 


-0*12 f B 


pv( 1 - W)B 


Diffusion of air 
from surface 


Convection 
toward surface 


(A2) 


since there is no net transport of air into the body in the steady 
state. From equations (Al) and (A2) , 


pvB = m(A + B) (A3) 

The equation for transfer of heat within the porous wall is 

(h^A + kB)^ = c^iCT - T c )m(A + B) (A4) 


Conduction into 
solid and coolant 


Rate heat absorbed 
by coolant 
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Equation (A4) may be written 

“ c pA< t - T =>* 


(A5) 


where k 
kjjA + kB 

A + B 
as 


is a mean value of the effective thermal conductivity 
Similar ly, equation (A2) may be written (by using eq. (A3)) 


pDnp — = mW (a6) 

dy 

— pD-i oB 

where pD-, 0 is a mean value of — =a_ . 

12 A + B 

Equations (A5) and (A6) have solutions (similar in nature) of the 

form 


T = T c + (Ty - 
which satisfies the condition T -» T c 


W =1 - (1 - 

A characteristic thickness 9,p may be defined as 



in which the temperature changes from its value at the surface T v to 
its value at the far interior T c . 

Similarly, the concentration W changes from W w at the surface 
to the value 1 in a layer of thickness 6^ (within the porous wall) 
where 


Ac p,iy 

T c )e * (A7) 

as y -* and 

J-y 

W„)e pE l£ (A8) 
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f 00 1 - ¥ a PD 12 

8 w = / itit dy = — 

v 0 w m 


(A10) 


B 


It Is seen that both thicknesses depend on the ratio — . It is desirable 

A 

to limit the relatively high temperature region within the wall to a 
narrow region near the surface; the specific heat Cp must be large 

and k must be small, that is, both k^, and k should be small. Con- 
siderations of the boundary-layer flow outside the wall also lead to the 
conclusion that the specific heat should be large and k should 

be s mall (so that the Prandtl number may be as large as possible). Thus, 
these requirements result in two desirable features: 

(a) Confinement of the high temperatures within the body to a 
narrow region near the surface and 

(b) Max imum convection of heat away from the stagnation point (as 
shown in the body of the report) 

The conditions at the surface as it is approached from within the 
wall are given as follows: 


From equations (A2) and (A3) 

-(pD^ ^ B = (1 - W w _)m(A + B) 
and from equation (A4) 

(kpA + KB)(^ = Cr,.l(T w - T c )i(A - 


(All) 


(A12) 


where the subscript w- denotes that the surface is approached through 
negative values of y. The rate of heat transfer from the stream 

bo un dary layer is fk — ] (A 4 - B) so that the boundary condition 

V <Sy/w 

(eq. (A12) ) becomes 

(k = Cp^(T w - T c )m (A13) 

Thus, there is a discontinuity in ^22 at the surface when k^. ^ k^. 
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Similarly, the condition of no net transfer of air to the body is 
given by 


"(p-^12 Tj ( a + B ) = U - W w )fvv w (A + B) (A14) 

\ 4y/w 

obtained when the surface is approached through positive values of y. 
Since 

Pv- V v- - * - Pw v w . 
equation (Al4) may be written as 



(l - W w )m 


(A15) 


It is noted that v is discontinuous since there is a discontinuity 


in the area available for convection given by the ratio 
concentration gradient 


B 


The 


dW 

dy 


A + B 

is also discontinuous for the same reason; 

dT 


this is expected by analogy with the temperature gradient Equa- 

ay 

tions (A13) and (A15) are seen to be identical with equations ( 3 ) and 
( 2 ) when evaluated at the wall. 
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Figure 2.- Boundary-layer mass balance. 
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Figure 3.- Effect of air-to-air injection on heat-transfer parameter. 
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Figure 4.- Concluded. 
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(b) Two-dimensional flow. 
Figure 11.- Concluded. 
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